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Abstract 

The Cauchy problem for the inelastic Boltzmann equation is studied for small data. 
Existence and uniqueness of mild and weak solutions is obtained for sufficiently small 
data that lies in the space of functions bounded by Maxwellians. The technique used 
to derive the result is the well known iteration process of Kaniel & Shinbrot. 

1 Introduction 

The theory developed by DiPerna & Lions in the 90's [13] on what is called Renormalized 
solutions has been a great success in finding existence theorems for the Boltzmann equation 
(BE): The Cauchy problem in |13j and the Boundary value problem [H]. The theory is 
strong and flexible and can be adapted to find solutions for different problems, for instance: 
The Vlasov-Poisson-Boltzmann system (VPB) [23], the treatment of the BE with infinite 
energy [2^, the relativistic Boltzmann equation and others. Indeed, a great deal of 
applications of this theory has been written in the last 18 years on BE related problems. 

The theory is based in the bounded entropy which is a feature of the elastic BE solu- 
tions. Unfortunately, it is not known how to obtain an a priori estimate that confirms such 
feature for the inelastic BE solutions, even in the case where a cut-off is imposed to the 
collision kernel. This simple fact creates a big upset for the theory in the inelastic case. 
It is important to say that in the one dimensional case, Benedetto &: Pulvirenti overcame 
this problem in [4J provided that the initial datum is essentially bounded and has compact 
support in the velocity space. The proof uses an iterative process introduced by J. M. Bony 
that relies strongly on the dimension. The technique allows to find a uniform control on 
the entropy, and hence, to prove existence and uniqueness provided the initial datum has 
the afore mentioned properties. 



It is clear that more understanding in the collision operator is needed to solve the in- 
elastic BE in its simplest form: The Cauchy problem. More complex problems, like the 
initial/boundary value problem or the VPB system, are still out of hand. This paper re- 
turns to the late 70's and presents an application for the inelastic Boltzmann problem of 
the technique introduced by Kaniel & Shinbrot in that time. Known as Kaniel & Shinbrot 
iterates [23], this technique was created by these authors to find existence and uniqueness 
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of solutions for the BE in sufficiently small time. The argument is beautiful and simple, 
however, the lack of a priori estimates makes difficult to use it for large time existence, 
particularly for hard potentials (see [1] for an interesting result). We refer the work of Ukai 
& Asano [29] for a successful extension from small to large time existence using Semigroup 
theory in the case of soft potentials. 

In the mid 80's, Illner &: Shinbrot |22j realized that it is possible to make a trade-off 
between the "size" of the initial datum and the size of the time interval where the existence 
of solutions can be proved. They modified Kaniel & Shinbrot argument to obtain global in 
time solutions in the elastic case. Indeed, Illner & Shinbrot [22] observed that if the initial 
datum was close to vacuum and had sufficient decay at infinity, the argument could be 
carried out after dominating the solution globally in time with appropriate estimates. See 
also Palczewski & Toscani [27 \ for an additional interesting application of Kaniel & Shin- 
brot iteration on the existence of the BE when the initial data is close to a local Maxwellian. 



Since the work of Illner & Shinbrot [22], an extensive study has been made in the cases when 
the initial datum is small or a perturbation of the vacuum. Regularity, asymptotic trend, 
stability and quantitative properties of solutions can be found in [2], [9], [l9], [20], and [28] . 
Recently, Glassey |18] was able to find such global in time estimates for the relativistic case 
using the space 

(l + |xx^|2)(i+^)/V <oo|. 

oo J 

The weight (1 + |x x more appropriate for the relativistic case, replaces the 

Maxwellian weight used originally by Kaniel & Shinbrot (see S.-Y. Ha [19] for different 
polynomial weights and L^-stability). The parameter < < 1 is related to the decay 
hypothesis imposed in the scattering cross section. It is important to observe that the 
advection coefficient in the relativistic case is bounded, this fact makes this case more han- 
dleable. However, the higher complexity of the collision laws introduces great difficulty in 
the calculations of the estimates. 



Finally, it is valuable to mention that for the inelastic homogeneous BE considerable work 
has been done, see [5], [6], [7J, [8j, [l2], [15], [E] and [25] for extensive studies. In the 
Maxwellian molecules case (the relative velocity in the collision frequency is approximated 
by the thermal speed) the reader may go to the work of Bobylev, Carrillo & Gamba [7] 
and investigate different aspects of this problem, for example, existence and uniqueness of 
solutions, self-similar solutions and moment equations. In the hard sphere case, one of the 
first complete studies of inelastic interactions was done by Gamba, Panferov & Villani in 
[17] . In this work the authors investigate the BE for diffusively excited granular media, i.e. 
they study the equation 

dtf-fiA^f = Q{f,f). 

Under precise conditions they find existence of a solution that becomes rapidly decaying 
and smooth for arbitrary small time. It is also proved that such solution has overpopulated 
high-velocity tails with respect to the Maxwellian distribution. In fact, they prove that the 
solution is controlled from below by 

i^exp(-a|e|3/2). 

For an extensive study of general hard spheres kernels in the inelastic homogeneous case, we 
refer the work of Mischler, Mouhot & Ricard [25]. In this paper the authors study stability. 
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existence and uniqueness for the inelastic homogeneous BE with general hard sphere col- 
lision rate. In addition, they investigate the long-time behavior of solutions, in particular, 
they show that solutions of the inelastic homogeneous BE problem collapse in the weak* 
measure sense to the Dirac mass after (and not before) infinite time. Therefore, a total loss 
of energy in the system occurs as i ^ oo. This collapse phenomenon was first observed in 
the late 90's by Benedetto & Caglioti [3] in a one-dimensional system of n-inelastic particles 
with constant restitution coefficient. It is interesting that this collapse does not occur in the 
Maxwellian molecules case when variable restitution coefficient with appropriate behavior 
for small relative velocities is used, see |7j for a detailed discussion. In [8j Bobylev, Gamba 
& Panferov develop a complete study of the high-energy asymptotic (energy tails) for dif- 
ferent inelastic regimes with forcing or heating terms. They complete this study by means 
of the Povzner-type inequalities technique. 

In general, all these previous studies of the space homogeneous problem are achieved using 
the weak formulation of the BE as the starting point, such approach is different from the one 
followed in this work for the space inhomogeneous case which uses the "mild" formulation 
of the BE. In addition, the difficulty introduced by the advection term ^ • V/ in the full 
BE used in the space inhomogeneous case will require a near vacuum assumption in our 
approach. Such assumption is not necessary for the homogeneous case that only requires 
an initial datum with mass and first two moments finite to find existence and uniqueness 
of solutions. 

This paper is organized as follows: Section 2 is devoted to present the problem and the 
notation. Section 3 presents the Kaniel &: Shinbrot iterates with small adaptations to the 
present case. Section 4 shows different estimates essential to find global solutions, and, 
section 5 presents the existence theorems and some conclusions. 

2 Inelastic Boltzmann Equation 

2.1 The Cauchy Problem 

Let us assume that we have a large space filled with particles that are considered as mass 
points. Assume that these particles are interacting with a specific law and that the particles 
are not influenced by external forces. A good model to represent such dynamical system is 
given by the equation 

|^+^.V/ = Q(/,/) in (0,+oo) xM"xM" (1) 

The function f{t, x), where {t, x, ^) G (0, cxd) x M" x R", represents for any fix time t and 
velocity ^ the distribution of the density of particles throughout space x, for the fix time t 
and the fix velocity ^. Thus, the physical meaning implies that / > 0. Equation ([1]) is known 
as the Boltzmann equation because it was derived by the first time by L. Boltzmann in 1872 
in his studies of dilute gases. Today, the BE is used to model not only dilute gases but 
statistical transport associated to dynamical systems that behave "like" a dilute gas. Some 
examples are large scale interactions of Galaxies, granular gases and chemical reaction gases. 

The Cauchy BE problem consist in finding a nonnegative function / such that it solves 
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the equation ([T]) and the initial condition 



/(O, e, x) = m, x) on {0} X X M" (2) 

for some nonnegative function /q. Note that the BE is a hnear transport equation equated 
to a term Q{f,f). This term is called collision operator and its purpose is to model the 
interaction or intramolecular interaction between particles. 

Some notation is needed before the collision operator is introduced. Write ^, ^* for the 
velocities of two particles just before they collide. The symbols are used for the veloc- 
ities of these particles just after the collision occurs. The law that relates pre-collision and 
post-collision velocities is given by 

1 + e 1 + e 
= ? ■ n)n and = -\ ^{u ■ n)n. (3) 

The variable u is used for the relative velocity between the particles 

and the unit vector n S S^~^ determines the impact direction, i.e. the unit vector that 
points from the ^-particle center to the ^*-particle center at the moment that the particles 
collide. The parameter e is called the restitution coefficient. The purpose of e is to describe 
the inelastic effect in the collision between particles. A good physical approximation is to 
take e as a function of the impact velocity, i.e. a function of |n • n|. Contrary to the elastic 
case (e=l), in the inelastic case (e < 1) the vector n does not bisects the angle between 
the pre-collision and post-collisional relative velocities. In this work it is assumed that the 
restitution coefficient is only a function of the impact velocity e = e{\u-n\). The properties 
of the map e : z — > e{z) will be stated carefully later on. 

In addition to equations ([3]), there are other ways to parameterize the interaction between 
two particles, see for example |25j for complete studies using the center of mass-relative 
velocity parameterization. However, the impact direction parameterization given by ([3]) 
is appropriate for the present work since it shows the dependence of the impact velocity 
explicitly. 

The notation with acute marks extends naturally to the pre-collision perspective using 
the symbols to denote the pre-collision velocities of a pair of particles. For example, 
from the pre-collision perspective, equations ([3]) read 

1 -I- 'e l + 'e 

i = 'i — Y~ ^'^ ■ ^* " — 2~ ^'^ ' ^^"^ 

where 'e = e{\'u ■ n\). Observe that in equations (HD the symbols represents the post- 
collision velocities of the particles. 

The Jacobian of the transformation ^ will be needed subsequently for the manipulation 
of the collision operator. It is not hard to realize that 

= e {\u ■ n\) + \u ■ n\ Cz {\u ■ n\) = 9z {\u ■ n\) , (5) 



J{\u-n\) 



4 



where 9 : z ^ z e{z). The symbols ez{z) and Oz{z) have been chosen to denote the deriva- 
tive of e and 9 respectively, avoiding the possible confusion with the post-variable symbols 

"e'(2;)" and ^^9'{zy\ In this way, the Jacobian is fully determined by the restitution coeffi- 
cient e and hence, the Jacobian is a function only of the impact velocity |M-n|. Consequently, 
we denote 'J = J{\'u ■ n\) and likewise for the symbols e' and J' . 

After agreed in the previous notation let us describe the explicit form of the collision oper- 
ator. For this purpose fix / and g two nonnegative functions, then Q{f,g) can be written 

as 

Q{f,g) = Q+{f,g)-Q4f,g) (6) 

where Qj^{f,g) is known as the positive or gain part of the collision operator and is given 
by the integral formula 

Q+{f,g)= [ I 7T7/('05(t)l«-"|dnde*. (7) 

The term Q-{f,g) is referred as the negative or loss part of the collision operator and is 
given by the integral formula 

Q.{f,g)=[ [ f{09mu-n\dndC.. (8) 

An important point to observe is that Q-{f,g) can be written in the simpler form 

Q.if,g) = f-Rig), 

where 

R{9)= I [ g{^*)\u-n\dnd^, = Cn [ g{^*)\u\d^* = Cn g * \^\. (9) 
Note that the constant C„ can be easily calculated as, 

Cn= [ \u- n\dn = 2 IS^-^I z{l - z'^f'-^'^l'^dz = , (10) 

7o n - 1 ' ' 

in particular for the three dimensional case one has that C3 = I-k. 

Remark : In the definition of the collision operator Q{f,g), the hard sphere kernel l-u • n| 
has been used. This is common in the modeling of inelastic granular flows since in these, 
the interaction between particles is in fact a physical collision, i.e. the particles are modeled 
as if they were hard spheres colliding. 

2.2 Restitution coefficient 

The assumptions on the restitution coefficient e are: 

(Al) z e{z) is absolute continuous from [0,+cx)) into (0, 1]. 

(A2) The mapping z —>■ 9{z) = ze{z) is strictly increasing. Thus, 9{z) is invertible. 
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(A3) Define for any /3 > the ratio 



exp ( -| 



exp (^-f 02 j 

Then, the fohowing integrability condition on -(/'^ must hold 

fl n-3 

(^;3(x) =2|5"-2| / V/3M(l--z ) ' tiz G L°°([0,+oo)). (11) 

JO 

The elastic case is obtained when e{z) = 1. Assumption (Al) is natural since physical 
restitution coefficients behave in a continuous fashion. In addition, this assumption allows 
computing a.e. the Jacobian of the transformation (jj]) using formula ([5]). 

Assumption (A2) implies, by equation ([5]), that the Jacobian is positive. Thus, the trans- 
formation dH) is invertible. Indeed, subtracting equations dU 

u' ■ n = —e{\u ■ n\) (u • n), (12) 

hence, it is straightforward to show that 

u' ■ n = — sgn(u • n)d{\u ■ n\) and u ■ n = — sgn(n' • n)9~''^{\u' ■ n\). (13) 

The explicit expressions of CiC^C'*) and are immediate from equations ([3]) with 

equations (fT3]) at hand. Also, observe that (A2) imposes an asymptotic pointwise trend on 
the restitution coefficient. Indeed, note that 

, , 9(z) 9(1) e(l) , 
e(z) = > = for z>l, 
z z z 

thus, the decay of the restitution coefficient is controlled by below with the function 1/z. 

Condition (A3) is an integrability condition imposed indirectly on the Jacobian and it 
is required for technical purposes. Although it may look strange, it is virtually satisfied in 
any practical case (for dimension n > 2). It is certainly satisfied trivially in the elastic case 
and the following examples of commonly used restitution coefficients. 

(1) e{z) = eo G (0, 1]. Constant restitution coefficient. 

(2) e{z) = jqr^^j where a > and 7 G (0, 1]. Monotonic decreasing restitution coefficient. 

(3) Continuous restitution coefficient that is elastic for small velocities and inelastic for 

large velocities. 

. . f 1 for Z < Zq r r. 

e(z) = < / N r ^ tor some zn > 0, 
[ eQ{z) for z > Zq 

where eo{z) is usually a smooth decreasing function that may go to zero as z ^ 00. 
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(4) Viscoelastic hard spheres (Brihiantov & Poschel [10]). In this case the dependence 
between the restitution coefficient and the impact velocity is given approximately by 
the equatioE0 

e{z)+az^/^e{zf/^ = 1, (14) 

where the parameter a > (being a = the elastic case) is a constant depending on 
the material viscosity. Let us verify that such restitution coefficient fulfill assumptions 
(A2) and {A3). Differentiating equation (fHI) with respect to z one obtains 

a/5 z^l^eizfl^ 

zez[z) - 



l + 3a/5zV5e(z)-2/5' 
therefore, 

, , , l + 2a/5 zi/5e(^)-2/5 



Z f P , 2 



In addition, for any /? > 0, 

M') = Jf^^^^-v[-^(.^■'-H4'))■ (15) 

From equation ()14p follows that e(2;) — > as z ^ oo, therefore 

— ^(2;)^ = (1 — e(z)^) ~ for large z, 

this suffices to prove that -i/'^ G ([0, +00)), which implies assumption (A3) for 
dimension n > 2. 

2.3 Function spaces and properties of the collision operator 

The natural spaces to seek the solution of the Cauchy Problem are those functions in 
:= L^(M" X M") that are bounded by Maxwelhans. 

M"''^ = {/ G : I/I < cexp(-a|xp)exp(-/?|^|2) for some c> O} , (16) 

where a and (3 are positive parameters. The notation M*^'^ is used when the function only 
has Maxwellian decay in the ^ variable. 

It is easy to check that the following one is a norm in these spaces 



a,l3 



:=inf {I/I <cexp(-a|x|2)exp(-/3|eP)}. (17) 



The notation adopted for the positive cone in M"'^ is 



|/ G M"'^ : / > 0} . (18) 



The subscript + is often used to denote the positive cone, thus, it will be used to denote 
the positive cone of any function space used in this work. 



^The restitution coefficient for viscoelastic hard spheres can be described as a power series of z^^^ as weU, 
i.e. e(z) = ak z'°^^% for details see [10] . 
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The Banach spaces LP{0,T; X) with X a Banach space and 1 < p < +00 will be used 
to handle the time variable. They are defined as the functions f{t) : [0, T] ^ X that satisfy 

\\f\\Lno,T-,x)--=[£\\fit)\fx) ' <+oo forp<+oo. 

This quantity is, of course, a norm for these spaces. The case p = +00 uses the same 
definition, with the integral changed for the essential supremum on [0,T]. As an example, 
take X = M"'^ to obtain the spaces L°°(0, T; M"'^). These are defined as 



L~(0,r;M"'/^) = |/:[0,r]^M"'^:ess sup « < 00 I . 

y *e[o,T] J 



(19) 



Similarly, the standard spaces C(0, T; X) and ^^'^(0, T; X) are used. Recall that C(0, T; X) 
comprises all continuous functions / : [0, T] — > X with the norm 



X ■ 



Also, recall that W^''^{{),T; X) is defined as those functions / G V'{Q,T]X) such that ft 
exist in the weak sense and belongs to L^(0,T;X). The norm used for this space is the 
obvious one 

ll/lliyi'l(0,T;X) = II/IIl1{0,T;X) + 11/* IIl1(0,T;X) • 

The collision operator defined by equations dS]),© and ([8]) enjoys some standard properties 
that are compiled here for the future use. 

Fix /, g G M"'^, then the gain and loss parts of the collision operator lie in L\_ and have 
the following properties: 

(Pla) /ign + ^Cjd^ — for any a, 6 G M. Conservation of mass and momentum. 

(Plb) /jjn f)\£,\^d£, < 0. Loss of energy and cooling effect. 
(P2) ||Q+(/,g)||^, = ||Q_(/,5)||i.. 
(P?,) ||Q-(/,5)ILi = IIQ-(5,/)ILi. 

The derivation of these properties is classical and follows after a change of variables {^, i^*} — > 
{'^/C*} ill the positive collision operator. The reader may go to [llj or [13j for the derivation 
in the elastic BE and to [15j for the derivation for the Enskog equation. 



3 Linear Problem 
3.1 Mild solution 

As it will be explained later in more detail, the existence technique of Kaniel & Shinbrot 
for the Cauchy problem of the BE consists in finding the solution / of the linear problem 

^ + e • V/ + Q_ (/, g) = /i with /(O) = /o (20) 

with g and h carefully chosen, so that, / approximates the actual BE solution. Problem 
(j20p is linear due to the bilinear character of the operator Q and the fact that g is fixed. 
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In order to solve problem (I20p it is common to introduce the "function along the tra- 
jectories" , 

f*{t,x,0 := f{t,Mi,x)), 
where ipt is the trajectory map of the transport operator T := dt + i - ^ 

Using this notation it is not hard to prove that if / is smooth, equation (j20p can be written 
equivalently as (recall that Q-{f,g) = f ■ R{g) where R{g) was given in equation ([9])). 

^(t) + f*R*{g){t) = h*{t) with /(O) = /o. (21) 

Equation ([2T]) is a good base to define the concept of solution because it does not demand 
the differentiability in the x variable for /, equation (j2ip does. Moreover, if / is smooth, 
equations (j20p and (j2ip are equivalent in the sense that / is a solution of the former if and 
only if is a solution of the later. In other words, equation (I2ip is a generalization of equation 
(El. 



Definition 3.1. Define mild solution in [0,T] as a function / G W^'^{0,T; L^), that solves 
a.e. equation in [0,T] x x W. 

Observe that for any fix t > 0, the Jacobian of the trajectory map is 1, 



dic,x) 



1. 



Thus, one has 



Lemma 3.1. Let f e LP{0;T;L^) with 1 < p < oo, then ||/*(t)||^i = 11/(011x1 for all 
t > 0. 

Note that the functions / and completely determine each other because the trajectory 
map V't is a bijection for all t > 0, therefore there is no ambiguity to express the subsequent 
results in terms of / or The choice between the former and the later will be made so 
that the expressions are maintained as simple as possible. 

Theorem 3.2. Let T > 0. Let /o G M^'^ , g G L°°(0, T; M°'^) and h > such that 

[\*{T)dT e L^{0,T;M^''^). (22) 
Jo 

Then ^21\) has a unique mild solution G L°°(0, T; M"^'^) given by the formula 

= /o exp J R'^{9){^)do^ + J /i*(r)exp ^- j R*{g){a)da^ dr. (23) 

Moreover, \\f{t)\\^i is absolutely continuous, and 

|ll/WllLi + IIQ-(/,5)(i)llLi = l|/^(i)ILi, (24) 
in particular f G C{0,T; L^). 
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Proof. Since g and h are nonnegative, we may take nonnegative smooth sequences {gn} and 
{hn} such that 

9n g and hn y h p.w. as n — > +00. 
Define the nonnegative sequence of smooth functions in the variable t 

f*{t) = fo exp ^ R*{gn){a)da^ + ^ /i#(r) exp ^ dr. (25) 

Clearly, 

/* f* P-w. as n ^ +00, 
where is the nonnegative function given by equation (|23p . Moreover, 

/o exp l\*{gn){a)da^ < /o G L'(0,r;Li) 

and, 

^ /i#(T)exp (^-^ "^"^-Z /i#(T)(iT G Li(0,r;Li). 

Therefore, the convergence — > is not only pointwise, is in fact in L^{0,T;L^). 

In addition observe that since the sequence {/n (i)} is smooth in time, by an elementary 
fact of ODE'S 



dt 



(t) + f*R*{gnm = h*{t) with /„(0) = fo. 



(26) 



Assumption (j22p implies that lv^{t) G L^(0,r;L^). Indeed, note that 



T 



h*{l 



dr 



h*{T)di 



<c«„ 



h*{T)di 



L°°(0,T;M"'/3) 



hence, 



/i#(t) ^ /i#(t) in LH0,r;L^). 
Similarly, after some easy computations one has the inequalities for < t < T, 

R*{gnm < R*{g) < Co (lbll^oo(o,T;AfO.^)) (1 + 1^1) 



and 



ft{t)< \ ||/olL« + 



a,f3 



h*{T)dl 



exp (-a|xp) exp (-/?|^|^) (27) 



L°°(0,T;M° 

for some nonnegative constants Cq and Ci depending in the quantities as stated. Thus, 

f*R*{gn)^ f*R*g in L\Q,T-L^). 
Equation ([26]) implies that 

dft 



dt 
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however, as it was said /if — > /* in L^(0, T; L^), thus it is concluded that 

df* 



c 



dt 



Therefore, and / belong to M^^'^(0, T;L}), i.e. / is a mild solution of the linear problem 
(pTD . Moreover, using dUD one concludes that f* G L~(0, T; M"'^). 



For uniqueness, it suffices to prove that when /o and h are zero, / is also zero. When 
/i = 0, problem (j2ip can be written as 



I (/#(t)exp Q\*ig)ia)da^^ = 
Assuming also that /o = 0, follows directly from this equation that = 0. 



Finally, it is a well known fact (see Theorem lA.ll in Appendix B), that for any real Banach 
space X one has the inclusion W^'^{0,T; X) C C{0,T; X). In particular, this inclusion is 
true for X = L^. Thus, integrating (j2ip in space and velocity and using the fact that 



df* 



dt 



d_ 

di 



f*{t) 



.. = d^ 11^(^)11-- 



the equation (p^ follows. This concludes the proof. 



3.2 Kaniel & Shinbrot Iteration 



□ 



Kaniel & Shinbrot method consists in building two sequences, {/„} and {un}, from the linear 
problem (I2ip . These sequences are built in such a way that one is monotone increasing while 
the other is monotone decreasing. The key point is to prove that both squeeze down on 
a mild solution of the Boltzmann Equation. The way to produce such sequences is the 
following: Assume that Iq, . . . ,ln-i,uo ■ ■ ■ ,Un-i are known, then /„(t) and Un{t) are the 
mild solutions in [0, T] of the linear problems 

dl* 

-^it) + Q#(/„,n„_i)(t) = Q#(Z„_i, ln-i){t) 
dt 



-{t) + Qt{unX~i){t) = Qf(n„„i,n„_i)(t) 



(28) 



with ln{^) = Wn(0) = /o- The construction begins with a pair of functions {Io,uq) satisfying 



n 

what Kaniel & Shinbrot called the beginning condition in [0,T], i.e. Mq G L°°{0,T; M_ 



and 



< lf{t) < lf{t) < uf{t) < u^{t) in < t < T. 



,# 



(29) 
#1 



Lemma 3.3. Let {Iq,uq) satisfy the beginning condition in [0, T]. Then, the sequences {In} 
and {un}, defined in f2S\). exist for all n and belong to L°°(0, T; M"'^). Moreover these 
sequences satisfy 

< ifit) < l*{t) <...< l*{t) < ...< u*{t) <...< uf{t) < u*{t) in < i < T, (30) 
consequently, {/« } and {n^} converge in M"'^ in [0, T]. 
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Proof. The beginning condition implements the first step in the induction. Next, assume 
if,-- - , if-i and u*,--- , all exist, belong to L°°(0, T; M'^'^) and satisfy for t G [0, T] 

0<l*{t)<---< lt,{t) < ut,{t) < ■■■< u*{t). (31) 
Since for any to functions g and / 

£/ < / if and only if g* < f* (32) 

it is concluded that 

0<Qt{lo,lo) < - - - < Qt{lk-i,lk-i) < Qtiuk~i,Uk-i) < ••• <Qf(no,no). (33) 
Using the estimate of Lemma 14.21 in the next section, 

f Q*{uo,uo){T)dT e C{0,T;M^^^), 
Jo 

therefore, using ([33]) 

/ Qtih-i,h-i){T)dT and / Q^(tife_i, UA._i)(r)(ir 

JO JO 

lie in C(0, T; M?'^). Moreover, it is clear that using in the sequence of inequalities ([HI]) 



< h-i < Uk-i < 



,# 



exp (-alx - t^P - . 



Therefore and Uk-i lie in L°°(0, T; M^'^). Applying Theorem 13.21 to obtain that l"^ 
and uf exist and lie in L°°(0, T; M"''^). 



Finally, observe that by §1^ and ([32]) 

< R*{lo) < ■ ■ ■< R*(lk-i) < R*{uk~i) < - - -< R*{uo). (34) 
Using (j33p and (j34p one concludes directly from (|23p that 

this completes the proof. □ 

Since and {un} are monotonic nonnegative sequences bounded by uq G L°°(0, T;L^), 
they actually must converge in L°°(0,T;L^) to some limits. Moreover, one has 

Lemma 3.4. Let /o G M"'^ anrf (/o,'Uo) satisfy the beginning condition in [0,T]. Denote 
the limits of {In} and {un} by l{t) and u{t). Then l{t) = u{t) for all t G [0, T]. Moreover 
this common limit is a mild solution of the Boltzmann equation 

^{t) + Qtif, fm = Q*{f, f){t) with /(O) = /o (35) 

with f* G L°°(0,r;M"'^). 
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Proof. The fact that l{t) = u{t) follows identically as in [23j, Lemma 5.2. It remains to 
show that the common limit is mild solution of (j35p . Thus, let f{t) := l(t) = u{t) and 
integrate equation (j28p to obtain, 



l*{t)+ Qt{ln{r)^Un-i{T))dT = h+ Q# (r) , (t) )dr. 



Send n —>■ +00 in this equation. Since the sequence {In} is increasing, it is possible to 
pass to the limit in the integrals to get an integral version of ([35]) . Thus, it is deduced that 
/# (z VFi'i(0,T;Li), and by means of estimate (gSD, f* G L°°(0, T; M^'^). □ 



Lemma 3.5. The mild solution of the Boltzmann equation found in the Lemma \3.4\ is 

unique in the space L°°{0,T; M^'^). 

Proof. Assume that and f^ are mild solutions of ([35]) that he in L°°(0, T; M"'^)+. Then 
for i = 1,2, 

f*{t)= [ Qt{hJi){r) - Q*{fiJ,i){T)dT, fi = fo. (36) 





Define j{t) := fi{t) — f2{t). After subtracting equations (j36p and taking the norm in 
space and velocity, 

Wmhi < [\\\Q+(.flJ){r)hr + ||Q+(j,/2)(r)|Li + ||Q_(/l,j)(T)||^, + \\Q^{j,f2){T)\\^,)dT 

Jo 

= 2 r(||Q_(i,/i)(r)||^. + ||Q_(j,/2)(T)||^0dT. (37) 
Jo 

The inequality in the right hand side of ([37|) follows after applying properties (P2) and (P3). 
However, 

IIQ-(i,/i)(T)||ii + ||Q-(i,/2)(r)ILi <C^max ff || (1 + |^|)j(r) \\^, . 

1=1,2 a,p 

Using last inequality in (|57|) 



^1 <C^max ff f\\{l + \^\)j{r) W^.dr. 

1=1,2 a,l3 Jo 



Using Gronwall's Lemma one concludes that j{t) = in [0, T]. We refer to [23J for additional 
remarks about uniqueness. □ 

4 Estimates on the Collision Operator 

Lemma 4.1. Let % and defined by the pre-collision formulas ^ then 
y exp [-a \x + T(g -' exp [-a \x + r(g -^g»)p^ dr < ^^^^ ^ exp (^-a . (38) 



Proof. Define the vector b by 



1 +'e 

b = S, = — i'u ■ n)n 
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therefore, 



1 +'e 

-b = — = — - — {'u ■ n)n. 



Thus, the fohowing equahty is obtained 



+ r(e -'Op + k + r(C = \x + Th\^ + \x-Th + 



TU\ 



x\^ + \x + Tu\^ + 2t'^ [\b\^ -h-u] . (39) 



Now, since u ■ n = —'e{'u ■ n) one obtains, 



\b\^ -b-u 



'u ■ n\ > 0. 



(40) 



As a result of (1391) and (SO 



J exp ^— a |a; + r(^ — '0|^^ exp 



exp(— a|x| ) / exp(—a\x + Tu\ ) dr. (41) 



Observe that. 



I , |2 I |2 I , , ,2 X ■ U 

\X + TU\ = \X\ — \X ■ U\ + \U\ I- — ; h T 



but — \x ■ > 0, then from the previous inequahty 



+ 00 

exp ( —a \x + Tn|^ ] dr < I exp ( —a \x + tu\^ ) dr 

— oo 



< / exp —a \u\ — — j — h T \ dr 



— oo 
oo 



/oo 
exp (— a|n|^r^) dr = — 
-oo 'a 



Using inequahty (|42|) in (|41|) one obtains 



(42) 



□ 



Lemma 4.2. For any 0<t <T and f* e L°°(0, T; M"'^) the following inequality holds 



f QtUJ){r) ^iT<fc„,^exp(-a|x|2-/3|e|2) /# 



L°°(0,T;M° 



(43) 



where 



The constant Cn only depends on the dimension. In other words, 



f Qt{f,f){r) dTeL^{0,T;M"'P). 
Jo 



(44) 
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Proof. Note that 

Qt{f,f){T,x,0 = Q+if,f){r,x + T^,0 = 

I I 7\jf*i^^^ + ^i^-'0,'0f*{r,x + T{C-'C*),'C*)\u-n\dndC. (45) 

hence, 



Q#(/,/)(r,x,0 



Using the fact that 



< 



# 



, , 777exp(-a|x + r(e-'0l^-/3ren 



exp (-a|x + r(e - /Sl'e^n |u • n\dnd^^. (46) 



exp (-/?re|2 - p) = exp ( - - P 



\u ■ n\ 



one gets from 

QtifJ)ir,x,0 <exp(-/3|eP) 



# 



L°°(0,T;M° 



/ exp(-/3|^,|2) 



1 / I -'e^ \ 
— exp -p——\'u ■ n|2 exp {-a\x + r(C -'^P) 

exp (-a|x + r(^ -'^*)p) |u • n\dnd^^,. (47) 
Integrating (j47p by r and using Lemma l4.1l it fohows that 



r Qt(/,/)(r,x,0 

JO 



vr 



dr < — — -exp (-a|x| - 



# 



L°°{0,T;M°^./3) 

1 - 'e2 „ 



'u • \u ■ n\dnd£_^. (48) 



Observe that due to condition (jlip 

1 / .l-'e2„ 



— exp -/?- 

gn-l fe 'J * 



As a result, it follows from (|48p that 



'n.nndn = </.^(|'n|)<||</.;3||j 



j[ \Qtif, /)(r, X, 1 dr < a^^l^ exp {-a\x\' - m') 



f# 



L°°(0,T;M° 



this proves the Lemma. 



□ 



5 Existence of Global Solution 

Theorem 5.1. Let T > and assume fo G M"'^ with 

1 



< 



where k^^p is given in estimate ^3^ . Then, the problem i35\) has a unique mild solution 
with f* G L°°(0,r;M"'^) nC(0,r;L^). 
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Proof. It suffices to prove that the beginning condition is satisfied globaUy. Let 

/^ = and = Cexp (-a|xp 

then 

l*{t) = foexp(^- 1^ R#(uo){T)dT^ and uf{t)=fo + I^Q*{uo,uo){T)dT, 



44^ 44^ 

clearly < Iq < 1^ 1^ ■ In addition, from the last expression and estimate (^3 



Note that 



u* 



<{t) < ( ll/oL,^ + ka,i3 u* ] exp (-a|x|2 - . 



= C, therefore it suffices to choose C such that 

a, 13 



\\foL,p + k^,(3C'' = C 
to satisfy the beginning condition. This is possible as long as 

this concludes the proof of the theorem. □ 

Theorem 5.2. Under the same conditions of Theorem \5.1\ the problem ^35\) has a unique 
weak solution in the space 



A=[f:f*eL^{0,T;M^^^)]. 



Proof. Note that the mild solution of problem (f35]l lies in L°°(0, T; M"'^). As a result, 
Q{f,f) lies in L^(0, T;L^). Therefore, one can use Theorem lA. 21 to prove that / is in fact 
a weak solution of ([35]) . 

The uniqueness follows because if / and g are weak solutions that lie in A then Q{f,f) 
and Q{g,g) lie in L^(0,T;L^) hence by Theorem IA.2I thev are mild solutions, thus invoking 
Theorem 15.11 it is concluded that f = g. □ 

Remarks : First, it is important to observe that the conservation of mass and loss of energy 
properties (-Pla) and (-Plb) were never used in order to prove Theorem 15.11 In this sense 
there is no hope to prove large data existence of solutions with this technique unless these 
properties are included in the argument. Indeed, in the theory of Renormalized solutions 
for elastic BE properties (-Pla) and (Plb) are essential (in addition to the bounded entropy 
property). 



Second, Theorem 15.11 shows that there is no clustering (accumulation of mass) in the solu- 
tion for small data as long as the restitution coefficient satisfies the conditions (Al ), (A2) 
and (A3) in section (2). Observe also that since G L°°(0, oo; M"'^) one has 

< / < Cexp {-a\x - t^P) exp (-/3|C|^) 
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for C > given in Theorem 15.11 Therefore, 



Urn f{t,x,^) = a.e. in R" x M". 

t^oo 



Thus, for sufficiently smah data the equihbrium is always the trivial one. Similarly as Illner 
& Shinbrot pointed out in [22j, it is possible to estimate the asymptotic decay for the spatial 
density 

pit,x) = I f{t,x,^)dC. 



Indeed, note that 



p{t, x)<C I exp {-a\x - t^p - < ^ [ exp (-a|zp) dz 



whence, 

p{t,x) = o(—\ as t ^ +00. 



Finally, observe that thanks to the properties (Pla) and (Plb) of the collision operator, the 
mild-weak solution / has the following usual properties for t > 0: 

(1) Conservation of mass and momentum 

/ / f{t){l+Odidx= [ [ foil +Od^dx and, 
Jr" Jr" Jr" Jr" 

(2) Dissipation of energy 

/ / fit) \i\''dUx < [ [ fo iCl^d^dx. 
In fact, the solution / has all moments bounded due to the Maxwellian decay in velocity. 

A Some facts for the mild and weak solutions for the Inho- 
mogeneous Boltzmann Problem 

This section includes some elementary results that are needed along the paper. Let us start 
with a well known result about the spaces Vl^^'^(0, T; X). Its proof can be found in |16j . 

Theorem A.l. Let X a Banach space, and let f G W^'^iOjT; X) for some I < p < +oo. 
Then 

(1) f G C(0,T;X) (after possibly being redefined on a set of measure zero), and 
(2) 

fit) = f{s) + / ftiT)dT forall 0<s<t< T. 

J s 
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In the same way that the mild solution was defined, one can define a different concept 
of solution for a transport problem, namely, the weak solution. Let T := dt + ^ ■ 'V the 
transport operator, O C xR" x R" an open set, and h G L^{0,T; Lj^^{^)). Then one can 
look for a function / satisfying the transport equation 

Tf = hm (0, T)xn 
/ = /o on {0} X n (49) 

in the following sense, 

Definition A.l. A function / G L^{0,T; Lj^^{Q)) is called weak solution of problem 
if for any V S D{{0,T) x n) 

- [ [ fT^l,= f f hi, and, 
JO Jo. 

(2) /(O, •) = /o a.e. in Vt. 



Jn Jo Jn 



The following is a classical result in linear transport that relates mild solutions and weak 
solutions 

Theorem A. 2. Take = M" x M" m problem (4^ ), and assume that /o G and h G 
L^{0, T;L^). Then, f is a weak solution if and only if f is a mild solution for this problem. 

Proof. Assume that / is a weak solution for problem (I49p . Thus, for any •0 G -D((0, T) x 



fTi;= hi;. 
Jn Jo Jn 

Let a G D{Q) and p G D{{0,T)). Take ij{t,x,(,) = p{t)a{x — t^,(,) in this equation and 
perform the change the change of variables (t, x, ^) {t,x + t(^, ^) to obtain 

/ f*it)p'mix,o= r [ h*{t)p{t)i;{x,i). 

Jn Jo Jn 

This works for any i;, hence 



- / f*p'{t) = [ h*{t)p{t)dt a.e. in Vi. 
Jo Jo 



Therefore, using the definition of weak derivative one concludes that 

^ = h#GL\0,T;L'). (50) 

Recall that / G L^{0,T; Lj^J, whence same property holds for f^. So for any compact set 
K C ^, one applies Theorem lA.ll to obtain that /* G C{0,T; Lj^J, and that for all t > 
the equation 

f*{t) = fo+ fh*{T)dT in (51) 
Jo 

holds. Since the right hand side of equation (I5ip belongs to -L^(0, T;L^), it must be that 
at each time t equation (j5ip holds in fact in L^, therefore, G L^(0,T;L^). As a result, 
/GH^i'nO,T;Li). 

For the converse use equation ()50p and proceed backwards using the same idea. □ 
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